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X-ray photo-electron diffraction is used to directly probe the intra-cell polar atomic distortion and 
tetragonality associated with ferroelectricity in ultrathin epitaxial PbTiOs films. Our measurements, 
combined with ab-initio calculations, unambiguously demonstrate non-centro-symmetry in films a 
few unit cells thick, imply that films as thin as 3 unit cells still preserve a ferroelectric polar distortion, 
and also show that there is no thick paraelectric dead layer at the surface. 
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I. INTRODUCTION 

Theoretical developments and novel experiments in 
the area of ferroelectrics have rapidly evolved over the 
last ten years, allowing further progress in the under- 
standing of this remarkable phenomenon. In particular, 
"nanoscale" ferroelectrics have attracted considerable at- 
tention Q, 0, 01 . The question of the existence of 
a critical thickness, in other words whether or not fer- 
roelectricity can be maintained at reduced dimensions, 
is amongst the most exciting topics of the field today, 
with very active experimental 0, 0, 0] and theoretical 

efforts [is mm- 

Probing ferroelectricity in thin films and nanostruc- 
tures is a difficult task, which requires advanced tech- 
niques. Among these, scanning probe characterization 
based on piezoelectric microscopy has allowed a ferroelec- 
tric ground state to be identified down to 40 A in thin 
Pb(Zro.2Tio.8)C>3 films 0], and x-ray studies on PbTiOs 
films suggested that 28 A films are ferroelectric 0. Di- 
electric and pyroelectric response measurements have al- 
lowed ferroelectricity to be identified in polymer films 
down to 10 A (two unit cells) [Fj. More recently, ul- 
trahigh vacuum scanning probe characterization based 
on electrostatic force microscopy was used to study fer- 
roelectricity in barium titanate nanowires with diame- 
ters as small as 10 nm 14] . On insulating substrates, 
lateral periodicity was observed via x-ray diffraction in 
thin PbTiOs films down to 12 A and attributed to al- 
ternately polarized domains 0,0]- In ah these studies, 
however, properties averaged over the complete ferroelec- 
tric structure were measured and no local information on 
the atomic displacements was obtained. 

In contrast, the photoemission based photoelectron 
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diffraction (XPD) used in our study presents two interest- 
ing characteristics: it is naturally surface sensitive, due 
to an electron escape depth of approximately 20 A (at the 
energy used here) ; and atomic displacements within the 
unit cell can be directly probed, allowing the non-centro- 
symmetric and tetragonal nature of the crystal lattice to 
be directly demonstrated. This turns out to be crucial 
for studying ultrathin films, and for discriminating the 
behavior of the surface from that of the body of the film. 

The paper is organized as follows. In Section II, we 
characterize the ferroelectric distortion of PbTiOs and 
contrast it to the natural atomic relaxation appearing at 
surface and interface. In Sec. Ill A, we discuss aspects 
of the measurement methods while, in Sec. Ill B we 
show the non-centro-symmetry of a 20 A thick film using 
XPD. The Sec. Ill C is dedicated to the tetragonality 
measurements as a function of film thickness with XPD 
using Pb as emitter, down to 4 A (one unit cell). We 
summarize and conclude in Sec. IV. 



II. FERROELECTRIC DISTORTION VERSUS 
SURFACE RELAXATION 

The studies were carried out on c-axis oriented per- 
ovskite PbTiO\3 ultrathin films epitaxially grown on 
conducting Nb-SrTi0 3 substrates. Above 490°C, bulk 
PbTiOs is a paraelectric insulator and a simple cubic 
perovskite structure with a lattice parameter of 3.96 A 
("para" -state, see Fig.^a)). In this structure, the Ti and 
O atoms are in perfectly centro-symmetric positions with 
respect to the surrounding Pb cage. At lower tempera- 
ture, the material becomes tetragonal and ferroelectric 
with a- and c-axis parameters of 3.90 A and 4.17 A, re- 
spectively 0, ^| , as illustrated in Fig. QJb) . The ferro- 
electric phase is characterized by a non-centro-symmetric 
structure where the O and Ti atoms are unequally shifted 
with respect to Pb. In a unit cell with the polar c-axis 
along the z-direction, O and Ti move either upwards 
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or downwards (with a larger O displacement) resulting 
respectively in a "down"- and "up" -polarized state, as 
drawn in Fig. Efb). 

In the surface region (five top unit cells) that is probed 
by the XPD technique, the evidence of such a polar 
atomic distortion could be the signature of a ferroelec- 
tric "up"- or "down" -state but may also arise from the 
natural atomic relaxation at the film surface and inter- 
face already present in the paraelectric phase. A proper 
interpretation of our data therefore requires independent 
quantification of both effects. To that end, a reference 
configuration ( "para-unrelaxed" ) is defined in Fig. Ola): 
it corresponds to the truncated bulk paraelectric struc- 
ture of PbTiC>3 with the in-plane lattice constant con- 
strained to that of SrTiC>3 (clsto = 3.90 A) and a conse- 
quent tetragonality co/asro — 1-03 [ij. Fig. 0b) ("ferro- 
unrelaxed") shows the atomic distortion of the "up" -state 
as determined for bulk tetragonal PbTiC>3 by Nelmes and 
Kuhs in Ref. [l5j with c = 4.17 A. In order to quantify 
the surface relaxation, density functional theory calcu- 
lations were performed within the local density ap- 
proximation (LDA) using the ABINIT 0] package. Two 
different supercells were considered : a thick PbTiC>3 
slab in vacuum (Fig. ETc)) and a SrTi03/(one unit cell) 
PbTi0 3 /vacuum stack|l3(Fig.Etd)). Insulating SrTi0 3 
was considered in our simulations since Nb doping is not 
presently affordable at the first-principles level |2Cj. To 
reproduce the substrate clamping effect, the in-plane lat- 
tice constant was fixed to the relaxed a-axis value of bulk 
SrTiC>3 |2l| . The atomic positions were then relaxed un- 
til the maximum residual atomic force was smaller than 
40 meV/A. Our calculations were restricted to (1 x 1) 
surface periodicity and did not allow for an eventual an- 
tiferrodistortive (AFD) e(2 x 2) reconstruction p i |23 | . 
The latter is not excluded but, as discussed later, was 
not evidenced on our films at room temperature. The 
distortions in the upper half of each supercell are re- 
ported in Fig. E[c) , (d) ("para-relaxed" state). For easy 
comparison with the experiment, because of the typical 
LDA underestimate of the lattice constant [2l| , the val- 
ues are given as a percentage of cq. The magnitudes of 
the ferroelectric and surface relaxation effects can now be 
compared. First, the cation-oxygen displacements due to 
ferroelectricity (11.6% - 8.3% of c , Fig.|l») are signifi- 
cantly larger than the displacements due to surface relax- 
ation/rumpling (3.4% - 1.4% of c , Fig. 0(c) and 3.3% - 
1.5% of c , Fig.|2jd)). Second, the mean layer displace- 
ment for the "up" -state (Fig. 0b)) and for the surface 
relaxation (Figs. (c) and (d)) are opposite. Third, the 
surface relaxation and rumpling effects are globally un- 
affected by the film thickness (Figs. (c) and (d)) and 
their amplitude decays very quickly in the interior of the 
film : they are already negligible two unit cells below 
the surface. This implies that the XPD measurements 
will be probing both the narrow relaxed surface region 
and a few unit cells below, essentially not affected by the 
surface relaxation. 



III. EXPERIMENTAL RESULTS AND 
DISCUSSION 

A. Experimental details 

The samples used in this study are epitaxial, c-axis ori- 
ented PbTiC>3 thin films grown on conducting (001) Nb- 
SrTiC>3 substrates using off- axis radio- frequency mag- 
netron sputtering |24l |25| . Topographic measurements 
using atomic force microscopy (AFM) showed that these 
films are essentially atomically smooth, with a root- 
mean-square roughness between 2 and 6 A over a 10 x 10 
/xm area. Room temperature x-ray diffraction measure- 
ments, for films with thickness > 28 A allowed us to pre- 
cisely determine the thickness and the c-axis parameter of 
the films, and to confirm their epitaxial "cube-on-cube" 
growth. 

After growth and characterization, the films were 
transferred ex-situ to a modified Vacuum Generators ES- 
CALAB Mk II photoelectron spectrometer. The XPD 
measurement system comprises a hemispherical electron 
energy analyzer with a three-channel detector, an x-ray 
photon source with two possible energies (Jw = 1253.6 eV 
and 1740 eV for MgKa and SiKa radiation, respectively), 
and a computer-controlled two-axis goniometer capable 
of rotating the photoelectron emission angle over the full 
hemisphere above the surface |26l |27| . 

The local geometry around a selected atom can be 
probed by performing an intensity versus emission-angle 
scan of a chosen photoemission line. Because of the chem- 
ical sensitivity of photoemission, a given atom type is 
then chosen by selecting one of its core levels. The outgo- 
ing photoemitted electrons exhibit a strongly anisotropic 
angular intensity distribution. This angular distribution 
is due to the interference of the directly emitted pho- 
toelectron wave with the scattered electron waves. The 
analysis of the interference (or diffraction) patterns is 
facilitated by the so-called "forward focusing" effect tak- 
ing place for photoelectron kinetic energies greater than 
w 0.5 keV. When considering a row of atoms, scattering 
at the first few atoms along this row focuses the elec- 
tron flux in the emitter-scatterer direction (for a review 
see Ref. [2^, |2j|). This enhancement of the intensity 
in the emitter-scatterer direction is schematically illus- 
trated by the green curve in Fig. 0[a) (right part) for 
the centro-symmetric "para" -state (continuous line) and 
the "up" -state (dotted line). The forward focusing ef- 
fect is further amplified for electron scattering by heavy 
atoms. In a semi-classical picture this can be understood 
as the focusing of the electron wave by the high number 
of protons in high atomic number atoms [2^. Note that, 
despite the forward focusing effect, the experimentally 
measured angles are sensitive to multiple interferences, 
refraction and possible anisotropic atom vibrations at the 
surface. In the present case of PbTiC>3, Pb scattering 
is highly dominant compared to the scattering by other 
elements [3(J. As a first step, in order to probe the non- 
centro-symmetry, O was chosen as emitter-atom (O Is 
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core level, E^j„ = 724.1 eV), since it has the largest dis- 
placement LlJ and has Pb scatterers as nearest neighbors 
(see Fig. El . However, the O contribution from the Nb- 
SrTiOs substrate becomes non-negligible for films thin- 
ner than the photoelectron inelastic mean free path, mak- 
ing the study of films thinner than 20 A more difficult. 
As a second step, in order to probe the tetragonality of 
the films, i. e., the c/a ratio of the Pb lattice (related 
to the polarization via the polarization-strain coupling as 
discussed below and in details in Ref . @) , Pb was chosen 
as emitter-atom (Pb 4f 7 / 2 core level, Efc in ~ 1115.5 eV), 
and Pb-Pb forward focusing directions were used. Since 
Pb atoms are absent from the substrate, this study can 
be done down to a monolayer of ferroelectric material. 



B. Non-centro-symmetric position of oxygen atoms 

First, considering oxygen as the emitter atom, fully au- 
tomated computer code for calculating electron diffrac- 
tion in atomic clusters (EDAC) via multiple scatter- 
ing [311, based on the muffin-tin potential approxima- 
tion |32| was used to calculate the XPD pattern. Fig.[3fa) 
shows four O Is core level emission (Ekin = 724.1 eV) in- 
terference patterns. One is the measurement made on a 
20 A thin film while the three others are multiple scatter- 
ing EDAC calculations of the "up"-, "down"- and "para"- 
state. Intensities are plotted in a stereographic projec- 
tion with the center corresponding to normal emission 
(polar angle 9 = 0°) and the outer border corresponding 
to grazing emission (8 = 70°). The strongest intensities 
(surrounded by red ellipses) correspond to the scatter- 
ing of O Is photoelectrons by Pb nearest neighbours (see 
Fig. ISl^b)). The white circle is a guide to the eye in- 
dicating the polar angle of maximum intensity for the 
measured interference pattern. It is evident that the po- 
lar angle position of this peak, which is directly linked 
to the O-Pb directions, is perfectly reproduced by the 
"up" -state calculation while the "para" - and the "down"- 
state simulations predict a different position. The "up"- 
state (down-shifted O position) corresponds to a smaller 
polar emission angle (9 up in Fig- EJb) ) appearing closer 
to normal emission (center of the interference pattern in 
Fig-EHa)). Such measurements have also been performed 
on films with thicknesses of ~ 500, 200, 100, 60, 44 and 
28 A, and all perfectly reflect the characteristics of the 
"up" -state p| . 

These conclusions, drawn from visual inspection of the 
interference patterns locally around the intensity maxi- 
mum (Fig.|3Ja)), are confirmed by a global matching ap- 
proach using a reliability (R)-factor to evaluate the qual- 
ity of the fit between the complete experimental inter- 
ference pattern data and theory (Fig. Bib)). The c-axis 
lattice constant value and the O and Ti shifts are the ad- 
justable structural parameters. A cut in the (100) plane, 
containing Pb and O atoms, is shown to facilitate the dis- 
cussion. In the calculation, O and Ti atoms are moved 
together and the dipole is continuously changed from the 



"down" -state to the "up" -state, crossing over the "para"- 
state. The best fit corresponds to the minimal R-factor 
value, which is reached when O and Ti atoms are shifted 
below the centro-symmetric position (parameters used 
for the "up" -state simulation in Fig. Ola)), with an R- 
factor value of « 0.34. In comparison, for the same c-axis 
parameter but the opposite O and Ti atom shifts (param- 
eters used for the "down"-state simulation in Fig.EIa)), 
the calculation gives a much higher R-factor of « 0.47. In 
between (zero O shift), in the centro-symmetric "para"- 
state, the R-factor is w 0.45 (parameters used for the 
"para" -state simulation in Fig. |^a)). 

It is important to note that surface relaxation and 
rumpling, neglected here, cannot weaken our conclusions; 
indeed they would give a picture resembling a "down"- 
state, the corresponding O-Pb atoms being shifted in the 
opposite direction than what is observed (see Fig. |2Ic)). 
Also, the possibility of a surface AFD reconstruction was 
explored without finding evidence for it in our room tem- 
perature experiments. 

This R-factor analysis therefore quantitatively con- 
firms the observations made in Fig.[3Ja), namely that the 
measured interference pattern is best simulated with the 
"up"- state. This demonstrates unambiguously that, for 
a film as thin as 20 A the O atoms have a non-centro- 
symmetric position in the Pb cage corresponding to a 
non-vanishing spontaneous polarization. 

Let us emphasize that piezoelectric AFM measure- 
ments performed after the XPD experiments on the 
thickest films (~ 500 A) confirmed the uniform "up" -state 
polarization while a uniform "down" -state polarization 
had been initially found for the same films just after 
growth. This confirms the monodomain character of the 
as grown sample and also indicates that the films are uni- 
formly switched from "down" to "up" -state by exposure 
to our conventional x-ray source, attesting for the switch- 
able character of the polarization. The details behind the 
switching are presently not known, but we believe that 
it occurs at the initial stage of the experiment while the 
measurement itself is essentially done in zero field. In 
fact, our results do not depend on the x-ray intensity, 
proving that the films are in equilibrium state during the 
measurements |34|. As discussed below, the agreement 
between the tetragonality deduced from x-ray diffraction 
and XPD also suggests that the measurements are per- 
formed in similar conditions. 



C. Tetragonality via lead emission 

In a second step, considering Pb as the emitter atom, 
XPD was used to determine the tetragonality. As demon- 
strated in Ref. |j, below 200 A the tetragonality de- 
creases as the film thickness decreases. This decrease 
is a consequence of the strong polarization-strain cou- 
pling in PbTiOa and a signature of a reduced polariza- 
tion in thin films. In Ref. [{j, this polarization reduc- 
tion was attributed to imperfect screening of the depo- 
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larizing field [l^. With XPD, using Pb as emitter, the 
tctragonality was measured down to the unit cell level 
as shown in Fig. 0J The absolute values of c/a, deduced 
from the forward focusing angles, are particulary large. 
This might reflect a strong enhancement of the polariza- 
tion in the probed surface region (of the order of 80% 
for c/a = 1.15, from the PTO polarization-strain cou- 
pling), even larger than in the theoretical prediction of 
Ref. [lC|. However, as previously stated, we are not nec- 
essarily measuring the precise atom-atom directions and 
the anomalously large forward focusing c/a might also 
be related to other effects (anisotropic atom vibrations 
at the surface, refraction and multiple scattering inter- 
ferences). Therefore a comparison with x-ray diffraction 
(XRD) 9] must be done at the relative level (Fig. |U left 
and right scale). 

To study the evolution of the tetragonality as a func- 
tion of the film thickness, the measured XPD values are 
compared to the c/a values obtained by XRD. The XPD 
measurement on Fig. 0] confirms the evolution of c/a ob- 
tained from x-ray measurements in Ref. [9j and agrees 
with the theoretical prediction (dashed curve) relying on 
the suppression of polarization due to imperfect screen- 
ing of the depolarizing field @- The similar thickness 
dependence for the XPD (very surface sensitive) and the 
x-ray measurements (average on the whole film) implies 
that the polarization evolves at the surface in the same 
way as at the interior of the film and that there is no 
thick paraelectric dead layer at the surface. In addition, 
the XPD tetragonality measurement shows a continuous 
decrease of tetragonality down to the thickness of one 
unit cell 0]. Two ribbons are drawn in Fig. 0] labeled 
with 1 and 2. They indicate the regions within which 
c/a values of 1.03 and 1.01 are crossed with respect to 
both c/a scales. For film thicknesses above two unit cells, 
the c/a values are larger than 1.03, the value expected at 
the bulk level for the paraelectric phase (resulting from 
the mechanical constraint imposed by the substrate, see 
also Fig.[!Ja)). This observation directly implies, via the 
polarization-strain coupling, that the films still have a fi- 
nite -although progressively reduced- spontaneous polar- 
ization. At thicknesses of one or two unit cells, as can be 



seen on Fig.0] c/a drops even more, reaching a value close 
to 1.01 for the one unit cell thick film 0|. This further 
decrease highlights that macroscopic elasticity no longer 
applies at such thicknesses where the interlayer atomic 
distances arc affected by surface relaxation and rumpling 
as shown by the ab-initio calculations (Fig. Gtd)). The 
measured tetragonality agrees with the computed value 
of 1.01 for the one unit-cell thick relaxed paraelectric film 
suggesting the absence of any additional ferroelectric dis- 
tortion at this thickness. 



IV. CONCLUSION 

This study thus directly demonstrates non-centro- 
symmetry, unambiguously a result of ferroelectricity in 
PbTiC>3 thin films down to 20 A. The measurements of 
the tetragonality, with a continuous decrease down to the 
bare substrate, show that even extremely thin films (3 
unit cells) have a c/a value larger than 1.03, attesting for 
the presence of a non- vanishing spontaneous polarization 
at this thickness scale. As the film thickness is reduced to 
a single unit cell, the experiments, together with calcu- 
lations, strongly suggest that both non-centro-symmetry 
and tetragonality are governed by surface effects, giving 
rise for our geometry to a polar relaxed structure but 
without switchablc ferroelectric distortion. 
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(a) "Para" Detector.; 




FIG. 1: (a) Schematic of the measurement setup. Core elec- 
trons have well-defined binding energies, and their photoemis- 
sion spectra exhibit characteristic emission lines. By selecting 
a particular emission line, photoelectrons from a given emitter 
can be chosen, thus probing the local real-space environment 
of the emitter. The interpretation is facilitated by forward 
focusing of electron flux along the emitter-scatterer direc- 
tion. The difference between two "forward focusing" peak 
positions is directly related to a modified emitter-scatterer 
direction, illustrated here by the "para"- toward "up"-state- 
O shift. Above 490° C, bulk PbTi0 3 has a cubic unit cell, 
with the O and Ti atoms centered in the Pb cage ("para"- 
state). (b) Below 490° C, bulk PbTi0 3 has a tetragonal unit 
cell: two equivalent ferroelectric configurations corresponding 
to two opposite polarizations, "up"-state and "down"-state. 
The displacements (in fractional units) from cubic phase sites 
are 0.111 and 0.037 for O and Ti, respectively |T^|. 
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FIG. 2: Schematic view of the atomic displacements resulting 
from ferroelectric distortion and from surface/interface relax- 
ation in the uppermost layers of a PbTiOa film epitaxially 
grown on SrTiOa. In bulk PbTiOa, epitaxial strain produces 
a tetragonality at zero polarization that can be estimated 
from the macroscopic elasticity theory as cq/clsto = 1.03. 
The ionic configuration resulting from the truncation of such 
a strained bulk paraelectric state (a) is considered as the ref- 
erence structure. Freezing the bulk ferroelectric distortion 
(as reported in Ref. [TElp into this reference structure re- 
sults in the "up" -state shown in (b). Additionally, the nat- 
ural ionic relaxation at the surface in the paraelectric state 
has been computed from first-principles (see text) both for a 
thick PbTiOs slab in vacuum (c) and a SrTi03/(one unit cell) 
PbTiC>3 /vacuum stack (d). Numbers in black correspond to 
the change of interlayer distances. Numbers in blue corre- 
sponds to the atomic rumpling in each layer (cation-oxygen 
distance). All the values are in % of Co, except those in italic 
that concern SrTiOs and are in % of asro- 
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FIG. 3: (a) Stereographic projection of experimental and the- 
oretical O Is emission line intensities on quarter hemispheres 
for four different cases: experimental data for a 20 A thin 
film and theoretical simulations taking into account "up"-, 
"para"- and "down"-state structures. Normal emission corre- 
sponds to the center of the plot and grazing emission (6 — 70°) 
to the outer border, (b) R-factor calculation to quantify the 
agreement between experiment and simulations for different 
structures, where a low R-factor corresponds to good agree- 
ment. A cut in the (lOO)-plane containing Pb and O atoms 
is shown to facilitate the discussion. A downward shift of Pb 
and O atoms optimizes the fit. 
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FIG. 4: Tetragonality (c/a ratio) as a function of film thick- 
ness for PbTiC>3 films (blue circles), for the Nb-SrTi03 sub- 
strate through PbTiOs films of one and two unit cells (green 
squares) and the Nb-SrTi03 substrate surface (red triangle). 
Black open circles are x-ray diffraction data and the black dot- 
ted curve is the result from model Hamiltonian calculations 
from Ref. Ribbons, labeled I and 2, indicate regions as 
discussed in the text. The left scale indicates the c/a val- 
ues extracted from the XPD experiment forward focusing di- 
rections using 8. The inset shows the crystallographic plane 
used to extract the c/a ratio. Pb-Pb scattering dominates 
compared to Pb-0 scattering allowing 8 to be determined via 
Pb-Pb forward focusing. 



